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Eigenphase shifts and eigentime delays near a resonance for a system of one discrete state and two 
continua are shown to be functionals of the Beutler-Fano formula using appropriate dimensionless 
energy units and line profile indices. Parameters responsible for the avoided crossing of eigenphase 
shifts and eigentime delays are identified. Similarly, parameters responsible for the eigentime delays 
due to a frame change are identified. With the help of new parameters, analogy with the spin 
model is pursued for the S matrix and time delay matrix Q. The S matrix is found to be put into 
exp[i(a -I- ba ■ h)]. The time delay matrix Q is shown to be given as Q = ^''''■(1 + Pa ■ cf + Pf ■ a) 
where the first term is the time delay due to resonance, the second term is the one due to avoided 
crossing interaction and the last term is the one due to a frame change. It is found that P^ + P^ 
equals unity. 
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I. INTRODUCTION 

Resonances observed in the energy-dependence of atomic and molecular processes often correspond to an autoion- 
ization or predissociation observed in photo-absorption, i.e., to the excitation of a quasi-bound state in a continuum 
region Their spectra are characterized by rich profiles providing a wealth of information on the atomic and 

molecular dissociation or ionization processes. Beutler's early observation of rare-gas levels auto-ionizing into a single 
continuum with angular-momentum conservation was reproduced by Fano through the formula |^ : 

^W = ^o^^^, (1) 
1 -t- 

where e is the reduced energy defined by 2(i? — £'o)/r; Eq and h/T represent the resonance energy and the mean 
life time of the quasi-bound state, respectively; ao represents the photo-fragmentation cross section to the continuum 
state; q is an index that characterizes the line profile (here —q is employed instead of the usual +q). Extension 
of Beutler-Fano formula to the dissociation into multi-channels was obtained by Combet-Farnoux The energy 
dependence of S matrix near an isolated resonance for such predissociation or autoionization into multi-channels was 
obtained long ago and well known By diagonalizing S matrix, eigenphase shifts 6i (i=l,2,...n) are obtained from 
its eigenvalues as 

S = Ue^'^U, (2) 



and utilized frequently as a tool for analyzing the resonances . Eigenphase shifts and corresponding eigenchannels 
are also extensively used in various forms in Multichannel Quantum Defect Theory which is regarded as one of the 
most general and powerful theories of resonance ^ . 

In contrast to photo-fragmentation cross sections and S matrix, the analytical formulas for eigenphase shifts in 
multichannel processes are not known. According to the numerical studies, eigenphase shifts in the multichannel 
system are known as showing complicated behaviors near a resonance due to the avoided crossing between curves 
of eigenphase shifts along the energy Such a phenomenon of the avoided crossing along the energy is not 

■ conspicuous in case of the formulas for photo-dissociation cross sections and S matrices. We report here the detailed 
studies on the behaviors of eigenphase shifts and also times delayed by a collision [pHlC|] for the system of one discrete 
state and two continua, which turn out to be described by functionals of the Beutler-Fano formula. Their behaviors 
in the system of one discrete state and one continuum are already studied . If one more continuum is added to that 
system, there enter new phenomena of avoided crossing between curves of eigenphase shifts and eigentimes delayed 
and of times delayed due to a frame change along the energy besides the resonance behavior. The avoided crossing 
interactions between eigenphase shifts along the energy always take place near a resonance if they are not excluded 
for symmetry reasons, as we will see in the later section. The addition of more continua may simply introduce more 
avoided crossings among the newly added eigenphase shifts. In this sense, the current system of one discrete state 
and two continua may serve as a prototype system for the study of the effects of avoid crossings between eigenphase 
shifts and times delayed due to the frame change on the resonance phenomena. 
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The import of the current work may be in the findings of the new parameters for the avoided crossing interaction, 
where S matrix can be put into a form of e*^ with H hermitian and the time delay matrix (or Ufe time matrix as 
called by Smith jl^]) can be put into a form of 1 + P ■ a where P is a polarization vector and a is the Pauli matrix. 
The analogy with the spin system under a magnetic field is exploited to find further properties of S matrix and time 
delay matrix. 

Section || describes the known properties of eigenphase shifts. Section 111 describes the form of S matrix near 
a resonance and the equation for eigenphase shifts. It also discusses some special cases and general characteristics 
of eigenphase shifts and eigentime delays. Sec. |^ obtains the formulas of eigenphase shifts as functionals of the 
Beutler-Fano formula for the system of one discrete state and two continua. Sec. ^ obtains the formulas of eigentime 
delays as functionals of Beutler-Fano formula. 



II. KNOWN PROPERTIES OF EIGENPHASE SHIFTS 



It is well known that , for an isolated resonance in single channel scattering, the energy dependence of a phase shift, 
S{E), is described by the following formula 

tanRi^)-^°] = ^(^ = 4 (3) 

where (5°, Eq and F denote the background phase shift, the energy and the width of the resonance, respectively. The 
second equality in Eq. (||) defines the dimensionless energy. It can be deduced from Eq. ^ that S{E) increases slowly 
except for the narrow energy range around Eg where its value undergoes a rapid change and is eventually increased 
by TT as energy varies from — oo to oo. Due to the property of tangent function, S{E) can be determined only up to 
UTT (n integer). It is usually taken to be 5° at E ^ — oo. Then, it becomes <5" + tt at i? ^ oo. In the absence of 
resonance, S is equal to J". The background phase shift S'^ is a slowly varying function of energy and usually taken 
to be constant as a good approximation. Thus, the increase of d{E) around Eq solely comes from resonance. Let us 
denote the difi^erence between S{E) and J" by SriE). With this notation, Eq. ^ becomes tani5r(-E) = tanJr(e) = 
— 1/e, and with the mentioned phase convention of setting S{E) to (5° at E = — oo, it can be transformed as follows 

Macek |^ generalized Eq. to multichannel processes as follows 

n 

2{E-Eo) = Y,^kCoi[Sl^S„,{E)], m = l,...,n, (5) 

k=l 

where SmiE) is the eigenphase shift defined in Eq. (|^), 5^1 is the background phase shifts in channel k, and F^ is the 
partial decay width to channel k. 

Hazi also showed that the eigenphase sum, which will be denoted as 6-s (E) , satisfies the "single-channel" formula 



k=l 



where F is the total width of the resonance, i.e. F = X]m=i ^m, and (5^ is the sum of the background eigenphase, i.e. 
(5^ = X]m=i ^m- The reduced energy e is the one defined in Eq. (|^) except that F now denotes the total width. 



HI. EQUATIONS FOR EIGENPHASE SHIFTS IN THE NEIGHBORHOOD OF RESONANCE 

The form of S matrix in the neighborhood of an isolated resonance in multichannel processes is well known and 
repeatedly derived in the past using various resonance theories . Though its form is basically same for all resonance 
theories, its detailed expressions may look differently from theory to theory. Here, Fano's configuration interaction 
theory is adopted for the treatment of resonance Q . In the configuration interaction theory, an isolated resonance 
in multichannel system is treated as regarding the system as composed of one discrete state (p and many continuum 
wave functions V'b^''^- The latter are taken to pre-diagonalize H. In other words, it is assumed that continuum 
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wave- functions are unable to interact directly with each other. They can, however, interact with each other indirectly 
via their interactions with the discrete state 4> if the following integrals 



{i^E^'\H\<^) , (7) 



are not zero. This means that H is no longer diagonalized in the combined space spanned by the basis set {(/i, i/;^'^"'''}. 

Eigenfunctions J'''^ of H in the space spanned by {(p , ip^^^'^ } can be obtained analytically as shown by Fano [^,1). 
By applying the incoming wave boundary conditions to those eigenfunctions as follows 



S matrix may be obtained B as 



j" 



E~Eo-i^Y.k\ykEV 



where S'j,^,, is the S matrix of the background scattering. Eq. (|^) is different from that of outgoing wave in that i 
is replaced by — i. The incoming wave instead of outgoing one is employed here as our interests are in the photo- 
dissociation processes. 

As stated in Ref. Ztt^^, |Vfc_Ep is the spectral width of the resonance peak and is denoted as F. According to 
Rcf. dHI], Eq. (H) may be rewritten as 

Srj - E U}^„(^^<"5,.,,. + ^c,.cy.)uj:,^. (10) 

In matrix notation, the above equation becomes 

5 = f/o [e-^'^" + icc) C/° = J7°AJ70, (11) 
where c is the column vector and its fc-th element is defined by 



Ffc is the partial resonance decay width to the background eigenchannel k. 

Let us consider obtaining eigenphase shifts 5 by the diagonalization of S matrix 

5* = Ue-^'^U. (13) 

If ^ = e^^"*" -f ice is diagonalized as 

A = Ve~^'^V, (14) 



(15) 



From Eqs. (|l|) and ([§,[/ = U^V . 

Let X denote the eigenvalue e"^"' of A, then the eigenvalue of yl = 
following secular equation 



ice can be obtained by solving the 



ic\ — X 



lC2Ci 

iczci 



lCiC2 



«CiC3 
JC2C3 



(16) 
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If we divide the first row of the above determinant by ici, the second row by ic2 and the fc-th row by ick in general, 
and if we introduce Xk defined by 



ICk 



the secular equation becomes 



2 

k 


X 


(or Xk 






Xi 


C2 


C3 ... 


Cl 


X2 


C3 ... 


Cl 


C2 


X3 . . . 



ICk 



(17) 



(18) 



The above equation for the system of n continua can be transformed into Sardi's form W3 



W{xk - Ck) + ^ 

k=l k=l 

By dividing Eq. ( p9| ) by its first term, one obtains 

n 

E 

fc=i 



(^k (*^m ^yn) 



Ck 



= 1. 



Cfe - a:*: 



(19) 



(20) 



By substitution of Eqs. (12) and (pTp, Eq. mt) becomes 



fc=i 



= E - En-i-. 



(21) 



Equating the real and imaginary parts, Eq. ( pl| ) yields two relations: one is Macek's formula and the other one is 
X]fe — T- This derivation of Macek's formula is simpler than the original one. The Sardi's form is given as a 
polynomial of order n and more convenient than Macek's formula in obtaining the solution. In particular, it shows 
explicitly that the number of roots of Macek's formula is the same as the number of continua, n. We will use it in 
the next section to obtain the solution for the system of one discrete state and two continua. Before doing that, let 
us comment on the general properties of eigenphase shifts obtainable from the Macek's formula. 

From Macek's formula, we could draw two properties of eigenphase shifts. Firstly, its differentiation with respect 
to energy yields 



dE 



E- 

.fe=i 



>o. 



(22) 



The above equation tells us that the first derivatives of eigenphase shifts with respect to energy near a resonance are 



£^o)/r| 



ooj are given as 



SI SI SI 5° 



up 



positive. Secondly, asymptotic values of eigenphase shifts (at \2{E 
to multiples of tt. 

Let us assume that (5j < (5§ < • • • < (5" and that they are all in the range of [0,7r], for simplicity, which can be always 
done without loosing generality by the Hazi's formula (||). Then, the above two properties on eigenphase shifts tell us 
that eigenphase shifts vary between two asymptotic values of their abscissas given by [^1,(52], [(^^ji^a], • • [<^°-i, <5,n]- 
If 5° becomes larger than tt, it is shifted by — tt and increases again toward Si- Since the sum of eigenphase shifts 
shows the resonance behavior and increases by tt around resonance, each eigenphase shift should similarly show up 
the resonance behavior between two consecutive asymptotes of the abscissas , (5°^]^], (i = l,...,n). When the 
consecutive eigenphase shifts are the same, let us say (5? = ^i+i, then one of two corresponding eigenphase shifts 
remains constant while the other one shows up the resonance behavior between two asymptotes [Sf , S'1^2] ■ 



IV. EIGENPHASE SHIFTS FOR THE SYSTEM OF ONE DISCRETE STATE AND TWO CONTINUA 



Let us now consider the system of one discrete state and two continua. For this system, Eq. (noh becomes 
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- ci){x2 - C2) + ci(x2 - C2) + C2(xi - ci)] = 0. 
i), we obtain 



Substituting Eq. (|T^) into (|2^), we obtain 

The roots of the quadratic equation ( p^ ) are obtained as 



■ 2\r -2i(5" I • 2\ I 2 2 n 



(23) 
(24) 



a; = e 




Ar f 

So) cos A?2 + — sin A?2 ± i < rir2 



E-Eo-i^ 

where T = Ti + r2, AF = Ti - r2, and A°2 ^ 5^ - S^. If we introdu( 
from the usual notation e in order to avoid the confusion with another 
later) as 



cos A5^2 



1° — (52. If we introduce a dimensionless energ; 
onfusion with another dimensionless energy 

E-Eo 




(25) 



(£;-£;o)sinA?2- — 

dimensionless energy (taken differently 
lensionless energy which will be defined 



the above equation becomes 



r/2 



e 



cos A"2 H — sin A j^2 i * 



— I 



r2 



^ . AO Ar 

Sm A"2 ^ COS A]^2 




(26) 



(27) 



where 5 is used instead of 5 in order to emphasize that the eigenphase shifts are now functions of instead of E. 

Analysis of the eigenphase shifts in Eq. (|2^) is done in the next subsection by introducing phase shifts due to the 
resonance and avoided crossing interaction. 



Eq. 



A. Phase shifts due to resonance and avoided crossing interaction 

|) can be decomposed into the product of two terms of unit modulus : 

Ar 



cos A52 



r 



■ sin A°2 T i 



4rir2 
r2 



• AO Ar ^0 

sm A"2 Y~ Ai2 



-I 1/2^ 



(28) 
(29) 



(the notations 5r and i5a without the bar are also used elsewhere if it is desired to express them in terms of E instead 
of Er). The phase shift 5r in the first term is the one due to the resonance as already mentioned in Eq. (^). In single 
channel scattering, it is the only term besides the background phase shift contributing to the eigenphase shifts near 
a resonance. In two channel scattering, there enters another phase shift 5a defined by the second term due to the 
indirect coupling of two continua via quasi-bound state near a resonance. 

The functional dependence of 5a{^r) on energy may be best seen by considering its cotangent like coiSj. — —Cr as 
in single channel scattering: 



Let us first rewrite Eq. 



cot 5a{tr) 



as follows 



Ur COS A°2 



AT sin A?, 



(e^ sinA?2 - cosA?2) 



4rir2 



r2 



1/2- 



(30) 



Tsin A?, 



cot 5a{^r) = — cot A"2 



2%/rTrJ 



AT 

r 



tan A?2) 



rsin Ajj 

"27rTrr 



Ur 



1fcotA?2) 



1/2- 



(31) 
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The term inside the bracket in the denominator of Eq. (|31|) can be transformed as follows 
rsinA° / Ar sinA°o / AF 



This suggests that we can introduce a new energy unit Ta'- 



cot A^2 = Lb - £;o - — cot Ay2 . (32) 



_ 2vrlr^ _ VP - Ar^ 

sinA02 " sinA02 ' ^ > 

For simplicity, let us also introduce a new parameter Ea'. 

Ar 

i;a = -Bo + — cotA;2, (34) 

(later it will be shown that Ea is the avoided crossing point energy and Ta is the strength of the avoided crossing 
interaction). 

With these parameters, Eq. ( |3^ ) becomes 2{E — Ea)/Ta which may be considered as a new dimensionless energy 
and will be denoted as e^: 

a 

By substituting Eqs. (||) and (|35|) into (|l]), we obtain 

cot ^.(e,) = - cot A;2 • ^^=, (36) 

where 5a{^a) is used instead of (Sa(er) in order to emphasize that it is now a function of and qa is the new parameter 
defined by 

-Ar 

gg EE ^ ^ -. (37) 

Vr2- Ar2 cosA?2 

Eq. (|3^ ) is a functional of the Beutler-Fano function defined by 

(e - q? 

/BF(e,<7)^^— ^, (38) 
1 + 

and can be rewritten as follows 



rnt A" \ - I cot A;2V7 BF(ea,ga) whcu < ,OQ^ 
COT Oa^eaj—< „ ^ , (^CSyj 

1^ -C0tA']'2V/BF(ea,ga) whcu > 

Eq. (|3^) or ( p9| ) tells us that the phase shift 5a due to the avoided crossing interaction takes its simplest form when 
it is parameterized in terms of and qa- Since we want Ta to be positive, it will be assumed that sin Aj'2 > 0, which 
can always be achieved by the appropriate choice of 1 and 2 for 51 and 52- With this convention in mind, Ta will be 
written as 2\/TiT2/\ sin A12I hereafter. 

With the above phase shifts, eigenphase shifts are obtained as 

25±{E)^5'i + 5l + 5r{E)±5a{E), (40) 

(for the choice of the sign, see ) . It is interesting to note that phase shifts of different origins are added up linearly 
to eigenphase shifts. 

Using Eq. (^), the above equation may be expressed in terms of the eigenphase sum 5^{E) 

25±{E)^5^{E)±5a{E). (41) 

The eigenphase sum 5^{E) consists of the background eigenphase sum 5^ and the phase shift 5r{E) due to the 
resonance. The phase shifts due to the resonance contribute positively from zero to tt as energy varies around a 
resonance. The phase shift is positive since particles are attracted and bound temporarily near a resonance. The 
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± sign in front of 5a{E) in Eq. (^) tells us that the phase shift due to the avoided crossing interaction in one 
eigenchannel increases at the expense of the phase shift of the other eigenchannel and vice versa. If it is increased 
in one eigenchannel due to attraction, it decreases in the other one due to repulsion. They are exactly canceled out 
and do not contribute to the eigenphase sum. The study of 5a as a function of shows that 5a has an extremum of 
cot~^(-y/(7^ + Icot A]'2) at ta — —^/<la minimum when go > 0, a maximum when qa < 0). From its definition, the 
limiting behavior of 5a at off-resonance is obtained as 



5ai(-a) 



^12 



A" 

^12 



when ea 
when e„ 



— oo, 
oo. 



(42) 



[note that < < 7r/2 according to the convention mentioned below Eq. (|39|)]. 5a increases around resonance 
by TT — 2Ai2- Thus the phase shift due to avoided crossing interaction varies most around resonance when — 
and least when = tt/2. It passes the middle, 7r/2, of two asymptotes of the abscissas in Eq. (B2) when = qa- 
Since the avoided crossing point energy Ea corresponds to zero of Ca [see the discussion below Eq. (|45D] , tells how 
the avoided crossing point energy is apart from the middle of two asymptotes of the abscissas. Such an apartness is 
a measure of the asymmetry of the curves of 5a{ea)- Fig. |l| shows the variations in the behaviors of 5a and 5± as qa 
varies. 

So far, we obtained eigenvalues of S matrix for two open channels near a resonance and decomposed them into 
contributions from the background, resonance, and avoided crossing interactions. Let us now consider obtaining the 
eigenvectors of S matrix which are usually called eigenchannels. S matrix is diagonalized in two steps. It is first 
transformed to A by which diagonahzes matrix, i.e., S — U'^ AU^ [S^ ~ C/^e^^'* U^). A is then diagonalized 
as in Eq. (Q) by V matrix which is composed of two eigenvectors f + and V- corresponding to J-), and 5- as 



After a lengthy derivation, eigenvectors are obtained as 

cos(^) 
sin(4) 



) 



-sin(%) 
cos(^) 



with 9a defined by 



cosf 



sin 9„ 



(43) 



(44) 



(45) 



00 
Thus, 



It is interesting to note that eigenvectors are independent of qa- They depend only on Eq. As varies from 
through zero to 00, 6a varies from zero through 7r/2 to tt and W4. varies from ^ q ^ through ^ |^ ^ to ^ ^ 

at Ca = 01 at E = Eq + -^cot Aj'2, two eigenphase shifts are avoided most. For this reason Ca = is considered 
as the avoided crossing point energy. (Another way to see that = is the avoided crossing point energy might be 
to remove the cause of the avoidance of two eigenphase shift curves and to let them cross through each other and to 
see whether the crossing point is ea 0. This was done and Eq = is confirmed to be the crossing point.) 
Let us parameterize U'^ matrix as 



[7° 



cos 6*° — sin 



sm ( 



cos ( 



(46) 



where 6*° may be considered as the angle by which the frame for background eigenchannels is rotated from the frame 
for background asymptotic channels. By calculating the multiplication of U — U'^V and by noting that U matrix is 
composed of eigenvectors m+ and U- of S matrix as 



U = {u+ U-) , 



the eigenvectors of S matrix are obtained as 



where 9'^ is defined by 



cos(-|-) 


, = 


-sm(^) 


_ sm(^) _ 




cos(^) 


0'a 


^Oa + 20° 





(47) 



(48) 



(49) 
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B. 5* matrix in terms of new parameters 



Let us now express S matrix in terms of new parameters introduced in the previous section. Reassembling by using 
its eigenvalues e~^^^+ and e'^"^^^ {6± — Ss ±Sa) and eigenvectors Eq. A becomes as follows 



_ -iSj: f COS Sa ~ i sin Sa COS 9a -ismSaSm9a \ ,t.„^ 
\^ —i sin (5(1 sin cos + i sin (5a cos 6*0 y ' 

From the relation S = U^AU^ with matrix defined in (Q), 5 matrix is obtained as 
^ _ / cos Sa - i sin 6a cos 0'a -i sin 6a sin 0'^ \ 



-i sin (5a sin 0^ cos ^a + i sin (5a cos 6*^ 

where 9'^ is defined in Eq. (^9|). Using Pauli matrices, it can be transformed into an invariant form: 

5 = e^**^ {cos 6a - i sin 6a {<Jz cos ^a + (^x sin ^a]} (52) 
= e^**^ (cos^d — isin(5a(T • ng/^) , (53) 

where the unit vector h^i is defined as 

fie;^ = icos6'^ +xsin0^. (54) 

Since 

cos — « sin(5a(7 • ne' —e *''<»'^"»i^ (55) 
it is simplified into the form of e'^ (iJ hermitian): 

(56) 

Since (t? • n^i^ — 1, a ■ hg/^ has two eigenvalues ±1. w+ and m_ in Eq. ( ^8| ) are their corresponding eigenvectors. 
In this case, it is well known that the vectors ingi^ equal the expectation values of the Pauli spin operator, 
(m±|(t|m±), and are called the polarization vectors that correspond to the eigenstates u±, respectively. We can obtain 
the change in the direction of the polarization axis as a function of energy by differentiating Eq. (^4|) with respect to 

^ = yxng,^-^. (57) 
dea dea 

Since tan 0a = — l/ca, 

dea dea 1 + Co 

The above equations also appear in the spin model for the adiabatic analysis of collisions |p^ . 

V. TIMES DELAYED BY COLLISION NEAR A RESONANCE 
A. Time delay matrix and eigentime delays 

For a single channel system, times delayed by collision, r, are obtained by the first derivatives of phase shifts: 

riE) = 2n'-^. (59) 



For a multichannel system, we may consider the time delay matrix Q defined by PJld 
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(60) 



The unitarity of 5* matrix ensures that the time delay matrix Q is hermitian. Let us now consider the time delay 
matrix for the system of one discrete state and two continua. Using Eq. (|5^) for 5" matrix expressed in terms of the 
Pauli matrices, the first derivative of 5* matrix with respect to energy is calculated as 

^ = -«-^5 + e "^"^ {-sm 6a- I cos Sa<T ■ ng'J — +e "-"^ I -tsm 5acr ■ j . (61) 

Now the time delay matrix Q is obtained as 

Since figi^ x y and y are perpendicular to he>^ , their linear combination is also perpendicular to figi^ . Let us introduce 
a vector 

nj^.^ = y X ng>^ cos Sa - ysinSa, (63) 
perpendicular to ngi and making an angle Sa with hg' x y and also define the following quantities 

= 2h—, 

rf = 2hs\n5a'^. (64) 

The last three quantities correspond to the time delay due to the resonance, the avoided crossing interaction, and the 
change of the direction of the eigenvectors figi of S matrix, or usually termed as the frame change, as a function of 
energy, respectively. Then, the time delay matrix Q may be written as 

Q ^]^{Tr + a ■ ngi^Ta + a ■ hj.jf). (65) 

Since observable properties of a particle with spin ^ can be represented as a function of its spin polarization P and 
the states of all the two-level (two-continua here) systems can be mapped on the states of orientation of a particle 
with spin i, Q may be rewritten in terms of the polarization vectors as 

Q = (l + Pa-S + Pf ■ a) , (66) 

where polarization vectors are defined as 

Pa^^ne-^, Pf = ^hj-,^. (67) 

The degrees of polarization \Pa\ and |P/| will be denoted simply as Pa and Pf. They are positive by definition (the 
term polarization may not be adequate for the time delay matrix. We will return to this point later). Let us introduce 
a new unit vector 

ht = figi I + fig, (68) 



Then Eqs. (|65|) and (|67D become 

Q^\(Tr + a- ht^rl + rf) - ^r, (l + • c?) , (69) 

where 



9 



Pt^^P^ + Pfht. (70) 
Now eigenvalues ol the time delay matrix Q or eigentime delays, which wiU be denoted as t±, are easily obtained as 

r± = i (r. ± ^T^ + r].) = Ir. (1 ± P,) . (71) 
Notice that Tq and Tf are added up incoherently to eigentime delays. Also because of tr(cri) — (z=x,y,z), 

tvQ = T++T^ =Tr, (72) 

that is, the eigentime delay sum has the same form as that of the eigentime delay in the single channel scattering [p^ . 

B. Explicit formulas of time delays 

The time delay due to the resonance, the first term of Eq. ([7l|), is obtained by differentiating tantf^ — — l/fr with 
respect to as 

..(..). 2.^^^. (73) 

It takes the Lorentzian form and is thus always positive or delayed as it is kept temporarily by the quasi-bound state. 
The form of the time delay due to avoidance, the second term of Eq. (|7l|), is rather complex: 

dSa [iAr(i? -i?o)sinA02 - (5)'cosA?; 



dE 



iE-Eor + {^)^] {rir2+ [(£;-£;o)sinA02-iArcosA02]'} 



By changing the independent variable from i? to e^, considering the reduced time delays Ta {—2h^^), and introducing 
the new parameters and q^- defined by 

2 r2-Ar2 



Eq. dTJ) is simplified as 



Sr^^cotAj^, (76) 



Tai^r) = —TrUr) — '^^ (77) 



The factor multiplying Tr{er) in the above equation can be transformed to the functional of the Beutler-Fano formula 
[the Lorentzian shape of 7v(er) may be regarded as a Beutler-Fano shape] 



Tri^r)—/ =, when e,. < g^, 

-Trier) I ^ =, when > g^, 



^^(^'■)\/ / BSe!';: fe^' whene, <g,, 

= Tr(er)5a(er)- (78) 

Interestingly, fa takes the Beutler-Fano formula in the energy scale of instead of Ca in contrast to the case of 5a, 
though the former is obtained as the derivative of the latter. Notice that the absolute value of ga is the polarization 
Pa considered in the previous section: 
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V /BF(er,gr) + 

ga changes its sign at — qr, so does the polarization vector Pa at the same energy. The study on Pa as a function 
of tr shows that Pa has a maximum of 1/(1 + q^)/ (1 + + r^) at — —^/qr and a minimum of at = qr 
(which corresponds to Ca = —1/qa where Sa has an extremum). Fig. |^ shows the variations in the behaviors of 
gai^r) and fa{er) as functions of at three different profile indices q^ with fixed (the value of q^. is restricted by 
|g,/cotA?2| > 1). 

Recalling that /BF(er, 9t) > 0, Eq. ( |79| ) tells us that the magnitude of Pa{er) is smaller than or equal to unity and 
has the effect of making the absolute magnitude of the time delay fa due to an avoidance smaller than or equal to 
that of the time delay fr due to a resonance: 

T.(e,) > |f,(e,)|, (80) 

which ensures that fr(er) i fa are positive. 

The above inequality may be obtained without knowing the explicit functional dependences of Tr and Ta as functions 
of energy. Notice that 

r.±r. = 2n^i%|M.2n^. (81) 
d-kj ah, 

Because of the inequality (p2|), the above equation should be larger than or equal to zero, which proves (|80|). This 
inequality is the manifestation of the physical fact that eigenphase shifts are increasing functions of energy in the 
neighborhood of a resonance as colliding particles are kept bound temporarily around a resonance (attractive forces 
make positive phase shifts). 

The above inequality restricts the magnitude of Pa to < Pq < 1. The magnitude of Pa causes the difference 
between two eigentimes delayed. Eq. ([7^) tells us that the degree of difference in two eigentime delays is governed 
by . Let us consider as a function of a; = Fi/F (0 < x < 1). Then since is symmetric with respect to a; = ^ 
and its values are zero at a; = and 1, infinity a,t x — ^, the magnitudes of Pa are 1 at a: = and at a; = ^. Since 
avoided crossing interaction is strongest at a: = ^ , the difference in two eigentimes delayed due to an avoided crossing 
interaction completely disappears when the avoided crossing interaction is strongest. 

Let us now consider obtaining the explicit formula of the times delayed due to a frame change. By substituting Eq. 
(p8|) and 



sin^a^^sinA;^^^, (82) 
into Eq. (|6^), the time delay Tf due to the frame change is obtained as 

Tf = Tr sin Aja W = ^rgf- (83) 

The absolute value of 5/ is the same as P/ : 

P/^lsinA^^lW^, (84) 




[actually g^ equals Pf in the convention mentioned below Eq. (|3^)]. Fig. || shows the variations in the behaviors of 
Pf{tr) and T/(er) as functions of energy at three different profile indices qr with A!j'2 fixed. 

As in the case of , Pf is smaller than or equal to unity. This can be proved by examining the behavior of the 
graph of Pf as a function of e^. P/ can be easily transformed as a function of by substituting 

F / AF \ 

cotA?2 , (85) 



for Ca into Eq. (p3|). It is then differentiated with respect to energy to yield the roots of the first derivative of ffi^r) 
at Cr = —1/qr and q^.. From this, it is found that Pf has the minimum of sin Aj'2^1 — (AF/F)^ at —l/q^. and the 
maximum of unity at qr- Its values become i/l — (AF/F)^ at ±00. This proves < P/ < 1. 
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Now let us consider the magnitude of the total polarization vector Pt = y-^a + ^/ • order to obtain this, let us 
fmther transform Eq. (|85|) using the definition (|33|) of Ta and AT/T = cot AJj/^t as follows 



|sinA;2 



cot A?^ 



Also after some manipulation, we obtain the following relation: 



cot^ A?2 



(86) 



cot 



a \2 



COt^ A?2 



^cotA?2^ 



l2 A?, 
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(87) 



Using Eq. ( |8^ ) to calculate + 1 and comparing the latter with the right hand side of Eq. ( |87| ) , we obtain 



1 = 



sm^ A?, 
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(88) 



where = (g,-/ cot Aj'j)^ — 1 is used. Substituting Eq. ( p8| ) into Eq. (| 
Beutler-Fano formula: 



From Eqs. 



V/BF(er,gr) + 

and ([zol), we obtain the surprising result: 



Pf can be put into a functional of 



(89) 



p! + Pf = 1- 



(90) 



This means that eigentime delays for the system of one discrete state and two continua are zero and Tr{er). Though 
time delays due to the avoided crossing interaction and frame change are asymmetric with respect to the resonance 
energy and therefore the energy of the longest lifetime does not match with the resonance energy, the energy of the 
longest overall eigentimes delayed is exactly matched with the resonance energy. 



VI. RESULTS AND DISCUSSION 



Let us summarize the results. In this paper, we considered eigenphase shifts and eigentime delays near a resonance 
for a system of one discrete state and two continua using Fano's configuration interaction theory. The eigenphase 
shifts are obtained as 



2^±(i?) = E^" 



Sr{E)±5aiE), 



where Sr{E) is the phase shift due to the modification of the scattering wave by the quasi-bound state and given by 
— arctan(l/er) and 6a{E) is the one due to the modification of the scattering wave by the other wave through the 
indirect interaction via the quasi-bound state and given as a functional of the Beutler-Fano formula: 



cot Sa{E) = - cot A?2 



f g - qa 
(1+ £2)1/2- 



In the above formula, energy is expressed in the unit of Ta — 2\/TiT2/\ sin A^al and shifted from the reso- 
nance energy Ea by ■^cotAj'a. The shape of the curve of 5a{E) is characterized by the line profile index 
qa = — 2Ar/(2-y/rir2 cos Aj'j)- The strength of the avoided crossing interaction felt by the eigenphase shifts is 
governed by the magnitude of Fa. No shift is expected when Fi equals F2, where the strength of the avoided crossing 
interaction is strongest. The maximum shift obtains when either Fi or F2 is zero. 

Eigenvectors v± oi A S — 5° A corresponding to eigenvalues 5±{E) are obtained as 



cos(%) 
sin(^) 



-sin(^) 
cos(%) 
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where 6a is defined by 



COS^a = ° , Sin^a = 



Eigenvectors arc independent of qa- They depend only on e^. Sinee qa stands for the apartness of the avoided crossing 
point energy from the middle of asymptotes of the abscissas of two eigenphase shifts, the above result means that 
the characteristics of avoided crossing interactions are independent of the asymptotes of the abscissas of eigenphase 
shifts. The corresponding eigenvectors of S matrix are obtained by replacing Oa with 0'^ = Oa + 2^o- 
With the new parameters, S matrix is found to be expressible as 

where ng^ = z cos O'^ + x sin 6'^ ■ In terms of Pauli matrices, the time delay matrix Q = ihS^ ^ is found to be consisted 
of three terms: 

Q = ^{Tr + ^ ■ no'Ja + ^ ■ ni'Jf)^ 

the one due to the resonance, the one due to the avoided crossing interaction, and the one due to the change of frame 
as a function of energy. Because of the last term, eigenfunetions of Q matrix are different from those of S matrix. 

The time delay due to the resonance takes a symmetric Lorentzian form and the time delay due the to avoided 
crossing takes a form of a functional of the Beutler-Fano formula: 



-Tr{E)—, i =, when < ^r, 

V (<!r--9x)^ 

Tr{E) — . , when e^. > Qr- 



The asymmetry of Ta as a function of energy is brought about by the nonzero value of qr which is proportional to 
the shift of the avoided crossing point energy from the resonance one. Thus the asymmetry of Ta is caused by the 
mismatch in the positions of the avoided crossing point and resonance energies. 
The time delay due to frame change takes the following form: 



c2 



Tf{E) = TriE)sinA'i2 



+ 1 



The above form may be understood from the fact that the changes of the frame spanned by the eigenvectors of S 
matrix is governed by 1 + while the time delay due to the resonance is governed by 1 + Therefore the ratio of Ta 
to Tr will be a function of (l + e^)/(H-e^). Tf{E) is also found to be transformed into a functional of the Beutler-Fano 
formula. 

In analogy with the spin ^ system, the time delay matrix Q may be expressed in terms of polarization vectors and 
the Pauli spin matrices as 

Q = ^Tr(l + Pa-a + Pfa^ , 

where polarization vectors are defined by 

Pa = —ne'^, Pf = 

Like the spin ^ system, it is found that the absolute values of Pa and Pf are restricted to < [Pal < 1 ^md < \Pf \ < 1. 
In the present ease a complete depolarization means that eigentimes delays are the same regardless of eigenchannels, 
while a complete polarization means that eigentime delays are and Tr{E) as a function of energy. Eigenvectors for 
eigentime delays due to an avoided crossing interaction and due to a frame change are orthogonal to each other and 

contribute to the total eigentime delays as y^r^ + r| = ^J~P^~+~Pj . It is found that P^ + Pj = 1 . This means that 

one among two total eigentime delays is zero while the other one is the same as the time delayed by the resonance state. 
Though time delays due to an avoided crossing interaction and a frame change are asymmetric with respect to the 
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resonance energy and therefore the energy of the longest Hfetime not matched with the resonance energy, the energy 
of the longest overall eigentimes delayed is exactly matched with the resonance energy. Though this is a surprising 
result, it should rather be so if we recall that all the partial life times obtained from the partial photo-dissociation 
cross sections as a function of energy are the same in the system of one discrete state and many continua. But the 
detailed study on their connection is beyond the scope of this paper. 

Present work reveals the dynamical parameters that govern the behaviors of eigenphase shifts and eigentime delays 
for the system of one discrete state and two continua. It may be applied to the systems of more than two continua 
by approximating that such systems are a cascade of independent two interacting continua. In that sense, current 
study may be considered as a prototype model for the "isolated" avoided crossing interaction. For the system of one 
discrete state and three or four continua, Eq. (^9|) becomes third and fourth order equation and its solution can be 
obtained by the Cartan and Ferrari's method. It may be highly desirable to do the similar studies on these systems. 
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FIG. 1. Sa{E) vs. ea and 5±{er) vs. tr are plotted for qa = 0, ±1, and ±3. Values for Si and 52 are 7r/3 and 7r/6, 
respectively. 
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Fig. 2 
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FIG. 2. 5a(er) and fa{er) vs. tr are plotted for three different profile indices qr = 0.6, 1, and 5 with A 
7I-/3. 
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Fig. 3 
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FIG. 3. gf{er) and ff{er) vs. er are plotted for three different profile indices Qt = 0.6, 1, and 5 with A?2 
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